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I calculate the QED coupling, a, directly in the MS scheme 
using an unsubtracted dispersion relation for the three light 
quarks, and perturbative QCD for charm and bottom quarks. 
Compact analytical expressions are presented, making this 
approach particularly suitable for electroweak fits. After 
a _1 (m T ) = 133.513 ± 0.026 is obtained in a first step, I per- 
form a 4-loop renormalization group evolution with 3-loop 
matching conditions to arrive at &' 1 (M Z ) = 127.934 ± 0.027 
for & a (Mz) = 0.120. The corresponding hadronic contri- 
bution to the on-shell coupling is Aa^ d (M z ) = 0.02779 ± 
0.00020. The error is mainly from m c , and from experimental 
uncertainties in e + e~ annihilation into unflavored and strange 
hadrons and r decay data. 

PACS numbers: 12.38. Bx, 12.20.Ds, 14.80. Bn. 



I. INTORDUCTION 

The increasing precision of experiments at CERN, 
SLAC, and the Tevatron calls for refined theoretical cal- 
culations with corresponding accuracy. For example, the 
precision of measurements of Mw, the Z width, and the 
weak mixing angle, s z , has reached or surpassed the per 
mille level. Most theoretical uncertainties are presently 
still negligible compared to the experimental ones. In 
contrast, the QED coupling constant (MS quantities will 
be marked by a caret), 



R{s) = 127rlmil (had )(s) = 
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escapes a precise theoretical computation from the fine 
structure constant, a, for \x > 2m 7r o when hadronic ef- 
fects must be included. On the other hand, knowledge 
of a(Mz) is indispensable for the extraction of the Higgs 
mass from precision data. In particular, with m t known 
independently from top quark production at the Teva- 
tron, Mw and s 2 z now serve as the most important Higgs 
probes, but are strongly correlated with a(Mz). Clearly, 
precise and reliable information on a{Mz) is needed. 

This has prompted a great deal of activity in the course 
of the past 4 years Jlr||. The traditional strategy is 
to exploit the analytic properties of QCD and obtain a 
subtracted dispersion relation (SDR) || for the on-shell 
coupling, a(Mz). By virtue of the optical theorem a 
weighted integral over the e + e~ cross section ratio, 



is obtained. The task is now reduced to construct the 
function R(s) for the entire energy regime s > 4m^ us- 
ing both theoretical and experimental information. This 
is a complex process, and the various papers 0-|5) differ 
by the methods by which experimental data are aver- 
aged, by the energy regimes in which perturbative QCD 
(PQCD) is trusted, by the treatment and parametriza- 
tion of resonances, etc. The situation is aggravated by 
the fact that the data sets in question often lack a thor- 
ough documentation, and much of the data taking was 
done at a time when per mille precisions were not an- 
ticipated. Consequently, an emancipation from old or 
imprecise data is indicated whenever possible. 

The key question is the energy domain in which the 
use of PQCD is adequate. PQCD and QCD driven data 
normalization for energies as low as a few GeV have 
been advocated by Martin and Zeppenfeld , but at the 
time encountered scepticism. A breakthrough came from 
three sides. First, theoretical calculations and mea- 
surements of the r lifetime have matured significantly 
over the last decade, with the extracted strong coupling 
constant a s , in perfect agreement with the value ex- 
tracted from the Z lineshape ||. Secondly, in Ref. || 
Hocker studies the invariant mass distribution in r decays 
following the spectral moments proposed by Le Diberder 
and Pich @. When varying his input data, he obtains 
impressively consistent results for the non-perturbative 
contributions to vector and axial- vector two-point corre- 
lation functions, which are moreover fairly small. Hav- 
ing passed this test, Davier and Hocker || finally fit to 
analogously defined spectral moments of R(s), with the 
conclusion that non-perturbative contributions to R(m T ) 
are negligibly small. Motivated by these developments, 
Kiihn and Steinhauser j|] present a state of the art anal- 
ysis of the SDR approach. 

In this letter I introduce a different method based on 
an unsubtracted dispersion relation (UDR). This is the 
natural framework for the computation of a(/i), defined 
at only one momentum transfer, q 2 = 0. By its very con- 
cept it is a purely perturbative quantity. Complications 
from non-perturbative physics arise only to the extent 
to which dispersion relations are used. This is hard to 
avoid for the three light flavors, but c and b quarks will 
turn out to be massive enough to be treated exclusively 
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within PQCD. The approach chosen in this paper also 
has an important practical advantage. Calculations of 
Higgs limits or \ 2 plots H require thousands of fits each 
with multiple function calls. A numerical (dispersion) 
integration within each call would be too expensive com- 
putationally, but within the approach introduced in this 
work, no numerical integration will be necessary. As a 
result a(Mz) can be self-consistently recalculated in each 
call, and the parametric uncertainty due to a s (which is 
a fit parameter) can be dropped. 



II. HEAVY QUARKS 

The polarization function in Eq. (Q) is defined through 
the current correlator, 

(w,-?V)n(9 2 ) = ' / rf 4 .x e ^<o|rj M (x) Jt/ (o)|o), 



where is the electromagnetic current. For a heavy 
quark it has been calculated up to 3-loop 0(aa 2 ) in 
Ref. m. The result for n(' l )(0) is expressed in terms 
of the quark pole mass. The coefficients grow rapidly 
and application to charm (bottom) quarks is impossible 
(questionable). However, the adverse coefficients are al- 
most entirely due to the employment of the pole mass, 
which is (due to quark confinement) not a well defined 
quantity It is therefore appropriate to reexpress 

nO)(0) in terms of the MS mass, m(/i), yielding (in agree- 
ment with Ref. [pi), 
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where Qh is the electric charge of the heavy quark, n q the 

2 

number of active flavors, and L = In Now all coeffi- 
cients are of order unity, indicating a reliable expansion. 
Moreover, all terms proportional to it 2 have cancelled. 
By setting fi = m(/x) the L terms can also be dropped. 
The remaining constant terms play the role of matching 
coefficients to be applied when the number of flavors in 
the effective theory is increased from n q — 1 to n q . This is 
familiar from the renormalization group evolution (RGE) 
of a s . Since n-loop matching must be supplemented with 
n + 1-loop RGE, inclusion of the O(ao 3 s ) beta function 
contribution is required and will be discussed later. 

Eq. (|3|) describes the contribution of a heavy quark in 
the external current. The n q — 1 light quarks appear- 
ing in internal loops must be treated as massless, since 
3-loop diagrams involving two massive quarks with dif- 
ferent masses have not been computed. It is indeed safe 
to neglect terms of O (a 2 rh 2 / to 2 ) , since in practice the 
heavy quark mass, rhh, is always sufficiently larger than 



all lighter quark masses, to;, and we will follow this ap- 
proximation throughout. Conversely, in 0(a 2 ) the heavy 
quark also appears as a loop insertion into a one-gluon 
exchange diagram (the "double bubble" diagram), and 
in the wave function renormalization of a light quark in 
the external current. The limit q 2 — > can only be per- 
formed when the heavy quark is decoupled, i.e., the a s 
definition for n a — 1 active flavors is used. This has been 
done in Ref. f|T< 
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Eqs. (|3) and (||) carry to 0(a 2 ) the decoupling of a heavy 
quark 14], as was first suggested by Marciano and Ros- 
ner |15] for the case of the top quark and generalized 
to 0(a s ) in Ref. fllq] . The same decoupling can also be 
applied to the MS definition of the weak mixing angle. 

With a s (rh c )/ir « 0.13 and the absence of non- 
perturbative effects, Eqs. (|^) and (Q) can be used reliably 
not only for b but also for c quarks. Complications with 
J/4" and T resonances are then completely avoided at 
the expense of the introduction of a stronger dependence 
on the quark masses. The numerical uncertainty due to 
rhb will turn out to be small, while rh c will introduce an 
error comparable to the one introduced through the J/ 1 ^ 
resonances in the SDR approach. 



III. LIGHT QUARKS 

I now turn to the three light quark flavors. Applying 
Cauchy's theorem to the contour in Fig. yields 

n(o) = i / — imn( s ) + l / -n( s ). (5) 

7T J S — It ZTTl J S 

Ami |s|=Mo 

The optical theorem applied to the first term, and the 
substitution s = ^e l9 to the second, brings the UDR 
into its final form, 



Mo 2tt 

n (3) (o) = T^ / —R(s) + ±[deh^(e). (6) 

12tt z J s — ie Zir J 

4raJ 

As in the SDR approach, the first integral can be evalu- 
ated using the measured function R(s) up to a scale fio 
where PQCD is trusted. Together with the second (called 
j(3) hereafter) this results for /io < Mj/^ in the 3-flavor 
definition a^ 3 '(fio)- Other values of /i are obtained using 
RGE, and other quark and lepton flavors are included at 
/.t = rfi(fi) using the matching description discussed be- 
fore. Special care is needed if /io > rh c , where convention- 
ally 4-flavor QCD is used. The clash with 3-flavor QED 
will generate some extra (non-decoupling) logarithms. 



2 



Im s 




- Re s 



Contour for an unsubtracted dispersion integra- 
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Indeed, following Refs. |g|JJ] I will use fio = 1.8 GeV and 
the result of Ref . S , as quoted in Ref . H , 
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R(s) 



(56.9 ± 1.1) x 1(T 4 . (7) 



The difference between this and the first integral in 
Eq. (H) (times Attci) can be neglected since, 
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The second integral in Eq. 
the help of Ref. pf, 
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|) can again be obtained with 

-4 C (3) + 2^%£>1 + 
Mo J 

C(3) + ?C(5)+F(4) 
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where in the 0(a s ) term I kept the small s quark mass 
effect (~ 2 x 10 -6 ). F(x) can be reconstructed from the 
absorbtive part of IT 3 )(s), i.e., R(s). Below threshold it 
can be well approximated as an expansion in x, despite 
the fact that rh 2 < Un, as can be shown by comparison 
with the exact result Jl7j (th e large quark mass expansion 
in O(a^) is also known fig] ). The coefficients in F{x) 
decrease even more rapidly than in R(s), 
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In principle, F(x) also applies to the tiny b quark contri- 
bution to IT 3 )(0), but without the non-decoupling loga- 
rithms in the first term. 



IV. RENORMALIZATION GROUP EVOLUTION 



The QED (3 function including QCD corrections reads, 



a 2 da , . 2 2^ n (0) 
dp z 
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where coefficients with a caret are scheme dependent. 
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and 62 are scheme independent and can be gleaned from 
the 3-loop (3 function for simple groups |L9]]. 8% is 
straightforwardly computed from Eq. ([}]), resulting in 



5 2 = N C C F Q 
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where for QCD we have N c = Ca = 3, Cp = 4/3, and 
Tp = 1/2. e 2 can be obtained from the first term with 
the substitution Q 2 „Cp — > Q 4 The coefficients in the first 



and the last term of 



Eq. © 



are familiar from 




On the other hand, the second term cannot be obtained 
from Refs. Q , since it cannot be disentangled from con- 
tributions with gluons in the external current. £3 has 
been obtained by Chetyrkin pc|] , and can also be recon- 
structed in the following way. Analytical continuation 
encodes the 4- loop order logarithms of 11(0) in 



where the r$ are the non-singlet coefficients (singlet con- 
tributions which are to be treated likewise are ignored for 
the moment). Denoting QCD j3 function coefficients by 
Pl 3 ^ , and the constant terms appearing in Eq. (^) in two 
and three-loop order by p 2 = CV (55/16 — 3£(3)) and p 3 , 
I find from a comparative analysis of leading logarithms 
in the SDR and UDR approaches, 
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S 2 _ R (3) 
— — T2 - P P2, 
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With the n from Refs. JT|,|l| I confirm Eq. @, and 
find, 
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where for QCD Tj.d abc d abc = 10/3. Eqs. @ and 
agree with Ref. p(J ■ Some of the terms in Eq. ( |l2| ) can 
also be checked with the four- loop QCD /3 function [B2| . 
There are delicate cancellations for n q — 4 and 5, for 
which 5 3 = —1.21 and 1.23, respectively (for 6 quarks 
5,3 = 5.79), and the next order might be larger with- 
out this being an indication of a breakdown of pertur- 
bation theory. By assuming no cancellations in 64, a 
conservative estimate of higher order RGE contributions 
is 1*4/5x1 < C%. With 
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& s(fJ>o) T , xr „(3) r 
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and the approximation, 

F = ln(l+X)<0(l), 



(13) 



Eq. (||) can be solved analytically, which is welcome for 
numerical implementations: 



/ft 



[(4 - <5 2 ci) (X + + 5x (c 2 - ca) ^- + S 2Cl Y - Sicj^-} . 

Effects from (3 2 ~ 2, and e 2 ~ 0.05, are well below 10~ 6 , 
and can be ignored. QED corrections to the QCD (3 
function induce a contribution which is formally of the 
same order as e 2 . It turns out to be negligible, as well. 



V. OTHER CONTRIBUTIONS 

To complete the calculation of a (/in), non-hadronic 
contributions have to be added to Eq. (£ 



-had}( ) = 
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The logarithms are the RGE effects of electrons and 
muons up to 0(a 3 ). For consistency, only the leading 
logarithms should be included in 3-loop order. Leptonic 
0(a 3 ) results are also available in the literature The 
third term is the corresponding matching effect analogous 
to the 13/12 term in Eq. (||). At this point I should stress 
that MS masses are used only as far as QCD is concerned; 
when quark effects (or leptons) get QED corrected, the 
mass is treated as on-shell. The last term is the bosonic 
contribution ( tne w± contribution to (3 is +7/4). 

0(a 2 ) correction are included for fermions, but not for 
bosons, because 2-loop electroweak calculations are gen- 
erally unavailable. 



VI. NON-PERTURBATIVE EFFECTS 

Thus far the discussion has been entirely within per- 
turbation theory. As stated in the introduction, non- 
perturbative contributions are expected to be small. In 
this section I will discuss the uncertainties associated 
with possible non-perturbative effects. 

The operator product expansion (OPE) [ ^5|j26| sup- 
plements perturbation theory with terms suppressed by 
powers of s. Dimension 2 terms can only arise from an 
expansion in the light quark masses. Therefore, D = 2 
operators can be treated perturbatively and the strange 
quark mass effect is already included in Eq. (|^) . Dynam- 
ical operators appear only at D = 4 and higher. For 
example, the strange quark and gluon condensates |26[ | 
are the dynamical operators of D = 4 and give rise to 
the contributions. 
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where the condensates are of order —m 2 K f 2 and Aq CD , 
respectively. Note, that these terms are suppressed by 
two powers of a s and therefore very small. They change 
a(m T ) by about —2 x 10~ 7 , an effect completely negli- 
gible. Effects from up and down quark condensates are 
suppressed by a further factor of m^/m 2 ^, and quartic 
mass terms are tiny, as well. 

As for the c and b quark contributions, Eq. (||) 
orginates from a heavy quark expansion and there is no 
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quark vacuum expectation value. There is, however, 
small contribution from the gluon condensate pflj2j] . 
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which is negligible for b quarks, and is between —3 and 
— 7 x 10 for c quarks, depending on the employed 
value for the condensate. This is still below other non- 
perturbative uncertainties discussed in the following. In 
conclusion, non-perturbative effects which can be de- 
scribed by local operators within the OPE are well under 
control. I will argue below, that this fact can also be used 
to limit possible other non-perturbative contributions. 

While the OPE takes phenomenologically into account 
a class of non-perturbative effects, its coefficient func- 
tions are still computed by perturbative means. There- 
fore, it cannot fully assess truly non-perturbative effects 
proportional to exp(— const/a s ). To discuss the valid- 
ity of the OPE requires an understanding of such effects. 
The leading correction to the perturbative treatment at 
short distances is believed to be associated with the one- 
instanton solution |2(||28| . In pure Q QCD the density for 
instantons of small size, p, is proportional to p8p9|] 



dp 
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p 2n+l VQ!{ 
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where the integer n > 2 can be fixed on dimensional 
grounds. Clearly, this density does not apply to large 
size instantons with p > Aqcd' which contribute unsup- 
pressed. For example, one expects a large size instan- 
ton contribution of order Aq CD to the gluon condensate 
(n = 2). We lack a full understanding of instanton ef- 
fects, but large size instanton contributions are described 
with the phcnomcnological parameters of the OPE. 

When the expression ([l7]) is integrated over p, one 
would encounter ultraviolet singularities for large enough 
n. This can be seen by changing the integration variable 
first to p — p^ 1 , and then to 



a a (p) 



(3) 
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Inserting the lowest order /3-function coefficient for pure 
QCD would result in ultraviolet divergences for n > 6 ]26| ] 
(see also Eq. ( p0| ) below). The interpretation of these 
divergences is that the integral is actually cut off at p ~ 
po, where po is the energy scale of the problem at hand. 
This limits the applicability of the OPE in two ways: (1) 
It introduces scale dependences in the matrix elements 
which are in conflict with the separation of short and 



x The inclusion of (almost) massless quarks rather extends 
the validity of the OPE 



large distance dynamics within the OPE. (2) Starting 
from some larger value of n there will be no suppression 
of higher order power terms. 

In order to be able to trust results derived within the 
OPE, one should therefore try to limit possible small 
instanton contributions. The small size instanton den- 
sity (^) can be integrated for small n, 



do: s _ 2 „ f ^ 

720 J a., P V a. 
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(19) 



where the A n are parameters of 0(1) (Ao — 1 by nor- 
malization). The a s distribution is centered around 
a s = f (1 ± -^=). For a s (M z ) = 0.120 this would corre- 
spond to instanton sizes p ~ 1.5 GeV -1 . However, the 
extra suppression factor p~ 2n ~ p 2n effectively shifts the 
distribution center to smaller values: 



= -(!■ 

3 n 



n , , 1 , 



(20) 



For large enough n the small instanton contribution will 
be dominated by the energy scale pq. One can thus read 
off the suppression factor for OPE breaking effects, 



A 



720a s 



0MA, 



(21) 



where A is a collective parameter again of 0(1), and 
a s = a s (p ). 

While one cannot estimate the parameters A n , one can 
try to put an order of magnitude bound on Ai using the 
phenomenological value of the two gluon condensate H , 



GG) « 0.04 GeV 4 . 



(22) 



I will assume that it is dominated by small rather than 
large instantons. This assumption will allow an order of 
magnitude bound on small instanton effects if one ignores 
the possibility of large cancellations. For n = 2 the inte- 
gral (|l9|) is contributed on average by instantons of scale 
pi ~ 0.9 GeV. Using that one can conclude, A 2 < 0.07. 

Even assuming that small instantons are (unlike the 
OPE power terms appearing in 1^) not suppressed by 
further powers of a s /ir, one would find from Eq. J2l| ) 
with A — 1 a suppression factor of 3 x 10~ 2 . From that I 
infer that non-perturbative contributions to or 1 should 
be < 0.006. 

By ascribing a discrepancy in semileptonic D decay 
data to OPE breaking instanton effects, it has been ar- 
gued J3(| that instanton contributions could be as large 
as about 10% in R{m T ). Contributions of this size are 
not excluded by the data, and would correspond to about 
100% of the PQCD corrections. Therefore, I will take 
50% of the PQCD correction to /( 3 ) as a (conservative) 
uncertainty introduced by possible OPE breaking effects 
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(which are expected to be insignificant away from the 
real axis). This yields an uncertainty of ±0.006 in a" 1 , 
incidentally identical to the bound obtained before. 

The SDR and UDR approaches are subject to the same 
size of non-perturbative effects, since we do not expect 
the relation between the on-shell and MS definitions of 
a(Mz) to be afflicted by low energy effects. Indeed, the 
authors in Ref. (3^] fitted different oscillating curves to 
the R{s) data around jj,o, and estimated the uncertainty 
to ±0.002 in a -1 . This is of the same order of magni- 
tude and smaller than the crude (and very conservative) 
estimate above. 



VII. NUMERICAL ANALYSIS 

In the numerical analysis I use a s (Mz) — 0.120 as 
a reference value. As mentioned earlier, no parametric 
error is included for a s , which is regarded as a fit param- 
eter. 

There is a variety of c |3^,|33| and b (3^]34|] running 
quark mass determinations. The uncertainties quoted by 
the various authors are of similar size, but being almost 
entirely theoretical, rather ad hoc. Therefore it does not 
seem justifiable to use a weighted average, which would 
yield0 m c (m c ) = 1.30±0.02 GeV, i.e. a very small error. 
Instead, I determine the averages and uncertainties from 
the spread of the results. This is a sclfconsistent treat- 
ment, as it only needs the usual assumption of normal 
error distribution^. The problem is then reduced to find- 
ing (posterior) information on a Gaussian distribution 
with (prior) unknown mean and variance, when given n 
data points (random drawings) y{. It can be shown pH ] 
that the marginal posterior distribution of the mean fol- 
lows a Student-t distribution, which is centered at the 
sample mean, 



1 " 

n ^ — ' 



(23) 



and has the standard deviation, 



1 n 
\ n(n - 3) ^ 



(24) 



Using this method I find, 



m c (m c ) = 1.31 ±0.07 GeV, 
rh b {rh b ) = 4.24 ±0.11 GeV, 



(25) 
(26) 



which introduce uncertainties of ±0.019, and ±0.002 in 
a^ 1 (Mz), respectively. These are added linearly since 
they are of similar origin. Notice, that the weighted 
and unweighted averages for rh c are (fortuitously) in very 
good agreement with each other. 

The 0(a 2 s ) term in Eq. is clearly dominated by the 

( 3) 

coefficient — /3q C(2)- I will therefore use the analogous 

coefficient in 0{af), -(/3[ 3) + 2/3^V 2 )C(2) » -19, as 
an estimate for the PQCD error in l( 3 >, corresponding 
to ±0.005 in a -1 . I note that this large coefficient is 
not problematic for PQCD, as terms of this type can be 
resummed to all orders. The total theoretical uncertainty 
is ±0.009. Eq. ([?]) adds an experimental error of ±0.015. 
A variation of a s within 0.120 ± 0.005 corresponds to 
a- 1 (M z ) = 127.934"°;°^, but this will not be included 
in the final error. 
The final result is 

a-^mr) = 133.513 ±0.015 ±0.009 ±0.019, 
a-^Mz) = 127.934 ±0.015 ±0.009 ±0.021, 



(27) 



where the errors are experimental, theoretical, and para- 
metric, respectively. (3! _1 (m T ) is to be compared with an 
earlier estimate ~ 133.29 |36| based on an O(a) calcula- 
tion and m t = 45 GeV. Using the 0(aa s ) conversion of 
Ref. |l(| (with an 0(aa 2 s ) improvement added), a(Mz) 
corresponds to 



Aa ( ^ d (M z ) = 0.02779 ± 0.00020. 



(28) 



Changing a s to 0.118 yields 0.02773 in perfect agreement 
with Aa^Mz) = 0.02774 ± 0.00017 from Ref. J§. Re- 
sults for other values of /i Z rh c can easily be obtained. 
A fit to all data using this approach yields for the Higgs 
mass, 



M, 



H 



compared to 



96j™ GeV, 



M H = 6912a GeV 



(29) 



(30) 



from a fit to the same data set ||, but using a(Mz) from 
Ref. |. 
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